Generalizing the differential algebra approach to 
input-output equations in structural identifiability 



Marisa Eisenberg 
Mathematical Biosciences Institute, 
Ohio State University 
and 

Departments of Epidemiology and Mathematics, 
University of Michigan, Ann Arbor 
marisaeOumich . edu 

Abstract 

Structural identifiability for parameter estimation addresses the question of whether it is possible 
to uniquely recover the model parameters assuming noise-free data, making it a necessary condition 
for successful parameter estimation for real, noisy data. One established approach to this question for 
nonlinear ordinary differential equation models is via differential algebra, which uses characteristic sets to 
generate a set of input-output equations which contain complete identifiability information for the model. 
This paper presents a generalization of this method, proving that identifiability may be determined using 
more general solution methods such as ad hoc substitution, Grobner bases, and differential Grobner 
bases, rather than via characteristic sets. This approach is used to examine the structural identifiability 
of several biological model systems using different solution methods (characteristic sets, Grobner bases, 
differential Grobner bases, and ad hoc substitution). It is shown that considering a range of approaches 
can allow for faster computations, which makes it possible to determine the identifiability of models 
which otherwise would be computationally infeasible. 

1 Introduction 

Identifiability analysis addresses the question of whether it is possible to uniquely recover the 
parameters from a given set of data. This problem can be broken into two broad (and sometimes 
overlapping) categories — practical or numerical identifiability identifiability, which incorporates 
the practical estimation issues of noise and bias, and structural identifiability, which considers a 
best-case scenario when the data are assumed to be known completely (i.e. smooth, noise-free 
and known for every time point). Structural identifiability is a necessary condition for parameter 
estimation with noisy data, and can yield information about how to reparameterize the model 
when it is unidentifiable. 

Many different approaches to structural identifiability analysis have been developed [IHl] ■ How- 
ever, the computational intensity of many methods makes applications beyond relatively simple 



models more challenging [S]. For linear models, identifiability can be determined globally via a 
transfer function approach and other linear algebra methods [H |6l [7] . One successful approach to 
identifiability for polynomial and rational function ordinary differential equation (ODE) models is 
via differential algebra [3 13 H] , which can be used to determine not only the overall identifiabil- 
ity of the model, but in the case of model unidentifiability, also uncovers identifiable parameter 
combinations and reparameterizations of the model in terms of these combinations j9] . In general 
these combinations are not unique and can be found in a range of different forms, although using 
Grobner bases one can find a 'simplest' set of combinations denoted the canonical set [9|. 

The algebraic approach is based on using characteristic sets [TU] (a method of solving/reducing 
differential polynomial systems) to generate a monic set of equations in terms of only the known 
or measured variables and the parameters, called the input-output equations, whose solutions are 
the set of all input-output pairs for the model. The coefficients of the input-output equations can 
then be used to test identifiability of the model [5j[Tll[T2]. One of the major limitations of this 
approach is that for more complex models it can become computationally intractable However, 
intuitively what makes the input-output equations informative about model identifiability is not 
that they are generated by characteristic sets, but rather that their solutions are the solution 
trajectories of the measured variables of the original system. It would therefore seem natural to 
extend this approach to more general methods of generating input-output equations. 

Indeed, an initial foray in this direction has been made by Meshkat et al. |13j, who show that 
Grobner bases provide an alternative to Ritt's pseudodivision in generating input-output equations. 
In this paper, we show that the input output equations contain complete identifiability information 
regardless of how they are generated, and moreover both Grobner bases (shown independently of 
[13]) and the differential Grobner bases of Mansfield [14] yield a set of input-output equations. Us- 
ing this more general approach allows more flexible calculations for testing structural identifiability, 
which can speed up computation and make tractable models which would otherwise be computa- 
tionally infeasible for the differential algebra approach. We then present examples and comparisons 
of the different solution methods, and show how choosing alternative solutions methods may allow 
the differential algebra approach to identifiability to be feasible for more complicated, nonlinear 
models. 



2 Identifiability 

We begin by introducing the identifiability framework used here. Let the model be given by 

x = /(x,t,u,p) 
y = 5(x,i,p) 

where x is a system of first order ordinary differential equations (ODEs), with t representing 
time, and u the experimental input function(s), if any. The model parameters are given by the Up 
dimensional vector p e M"p (the complex numbers C may also sometimes be considered, depending 
on the model). We will occasionally refer to individual parameters within p asp (without bold-face), 
and refer to an arbitrary point in parameter space (R"p) as p*. The measured data/output (s) are 
given by y, which represents the the riy-dimensional vector of output (s) without any measurement 
error. We also let Xq represent the vector of initial conditions for x(t). As in [HITT], we assume that 
/ and g are rational polynomial functions of their arguments, and that u, x, and y are arbitrarily 
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differentiable. We also assume that constraints reflecting known relationships among parameters, 
variables, inputs, and outputs are already included in the model equations, as these are known to 
affect identifiability properties [31 [71 [9] . Alternatively, the equality constraints among parameters 
or variables can also be appended to the model equations without affecting the results given here. 

Structural identifiability analysis explores the question: given an input u, model x — /(x, i, u, p) 
and experimental output y, is it possible to uniquely identify the parameters p, assuming "per- 
fect" noiseless data? Mathematically, this can be thought of in terms of injectivity of the map 
$ : p — >■ y given by viewing the model output y as a function of the parameters p [3 [9] . We note 
that because there may be some 'special' or degenerate parameter values or initial conditions for 
which an otherwise identifiable model is unidentifiable (e.g. if all initial conditions are zero and 
there is no input to the model), structural identifiability is often defined for almost all parameter 
values and initial conditions [3 [HI [H] . 

Definition 2.1. For a given ODE model x = /(x,t, u, p) and output y, an individual parameter 
•p is uniquely ( or globally) structurally identifiable if for almost every value p* and almost all initial 
conditions, the equation y(x, i, p*) = y(x, <,p) implies p = p* . A parameter p is said to be non- 
uniquely (or locally) structurally identifiable if for almost any p* and almost all initial conditions, 
the equation y{x,t,p*) = y(x, t,p) implies that p has a finite number of solutions. 

Definition 2.2. Similarly, a model x = /(x, t, u, p) is said to be uniquely (respectively non- 
uniquely) structurally identifiable for a given choice of output y if every parameter is uniquely 
(respectively non-uniquely) structurally identifiable, i.e. the equation y(x, t,p*) = y(x, t,p) has 
only one solution, p = p* (respectively finitely many solutions). Equivalently, a model is uniquely 
structurally identifiable for a given output if and only if the map $ is injective almost everywhere, 
i.e. if there exists a unique set of parameter values p* which yields a given trajectory y(x, i,p*) 
almost everywhere. 

The equivalence classes generated by $ are precisely the sets of parameter values yielding the 
same output, so that if the fibers of $ contain finitely many elements, the model is locally (non- 
uniquely) identifiable, and if the fibers of $ contain infinitely many elements, the model is termed 
unidentifiable. 

3 Algebra and differential algebra background 

In this section, we present a very brief overview of the differential and computational algebra con- 
cepts needed for this paper. For full details on the fundamentals of differential and computational 
algebra methods, the reader is referred to [TUl [HI [TS] . Let R be a ring in the usual algebraic sense, 
and x a set of indeterminates. For our applications, R will represent the field of coefficients for 
an ODE model, so that we will typically consider R to be either the real or complex numbers, 
depending on the model. We use the usual notation R[x] to represent the polynomial ring in x, 
and R(x) to represent the field of rational functions of x with coefficients in R. 

A differential ring is simply a ring in the usual algebraic sense, together with a differentiation 
operation which obeys the usual linear and product rule properties for derivatives. For ODE 
models, we typically extend R[x] to form a differential ring with derivatives in time t, denoted 
R{x}, by adding an additional derivative operation, in this case the usual polynomial derivative 
where we take derivatives from the ring of constants R to be zero. Elements of R{x} can be 
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thought of as elements of R[x, x', x", . ■ • ], where x' represents the set of derivatives of elements of 
X with respect to t (where x is our set of variables) . We note that for convenience we will often 
view a particular differential polynomial in x as an element of R[x, x', . . . , x''"''], where is n is the 
highest derivative of x appearing in the polynomial. 

Typically when working with differential polynomial rings, a ranking on the variables is chosen, 
in our case a ranking of the form u<u<u<---<y<y<y<... <x<x<x<-- - given in 
This makes allows one to determine leading terms, make polynomials monic, etc. The leader 
of a differential polynomial is defined as the highest ranking derivative of that polynomial (which 
can be a derivative of order 0). Choosing a ranking allows us allows us to fix the coefficients of 
the input-output equations uniquely, by dividing by the coefficient of the leading term to make the 
polynomials monic [U |9] . 

Let 5 be a set of differential polynomials in R{x}. The set of all polynomials that can be 
formed from elements of S by addition, multiplication by elements of R{x}, and differentiation is 
called a differential ideal generated by S, which we write as {S}. For a given set of polynomials 
S (or differential polynomials, where we simply view the derivatives of variables as additional 
indeterminates), the variety V{S) is defined in the usual way as the set of points for which all 
polynomials in S are zero. A differential ideal / is called prime if a6 G / implies that either a £ I 
or 6 G / and is called perfect if a'^ e / implies a G / (i.e. a perfect ideal coincides with its radical). 

There are several methods for manipulating systems of polynomials and differential polynomi- 
als, including the Grobner basis and characteristic set methods discussed here, as well as methods 
of resolvents, among others [HI [16]. The usual method used to generate the input-output equa- 
tions in identifiability for ODE models is the characteristic set [lOj . A characteristic set of a set of 
polynomials is defined to be a chain of minimal rank in the differential ring, where chains of poly- 
nomials are formed by using pseudoreduction |f 0] to reduce the rank of the polynomials compared 
to one another, until a minimal, autoreduced set is reached. For details on characteristic sets and 
their uses in identifiability, see [U |8[ [HI [17] . 

Grobner bases are one of the most common tools in computational algebra, and hence are a 
natural generalization of the characteristic set approach for generating input-output equations, as 
numerous fast methods for calculating Grobner bases have been developed (e.g. Faugere algorithm 
[18) . Grobner walk methods [H]). A Grobner basis of an ideal / is a generating set for that ideal 
such that the remainder of any element of ring yields zero if and only if that element is an element 
of /. Form more information on Grobner bases, the reader is referred to |15]. When applying 
Grobner bases to differential polynomials in x, we typically work in R[x,x', . . . ,x*-"-'] (where n is 
the highest order derivative appearing in our set of polynomials), where we treat derivatives of x 
as new indeterminates. 

ft would be natural to extend Grobner basis theory to the differential case, and indeed there 
have been two major formulations of differential Grobner bases [14, 16, 20 . Part of the difficulty 
in extending the Grobner bases to the differential case is in incorporating the differential structure 
of the ring, which lends itself to psuedoreduction rather than conventional reduction as is done in 
algebraic Grobner bases [Ml [16j . Mansfield Grobner bases [14] surmount this difficulty by devel- 
oping a pseudoreduction formulation of differential Grobner bases, and so we use this formulation 
here. A Mansfield differential Grobner basis of a differential ideal / is a generating set of / such 
that full pseudoreduction of any element of / yields zero. Details and comparisons of Mansfield 
and algebraic Grobner bases can be found in Mansfield's thesis [T4] . 
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4 Generalized differential algebra approach to identifiability 



Next, we introduce the basic differential algebra method, and give a few results which generalize 
the usual approach. As discussed above, the larger goal of the differential algebra method is to 
use algebraic manipulation to yield a monic set of Uy equations only in the input variables u, 
output variables y and parameters p, denoted the input-output equations [8], for which solutions 
are precisely the allowed input-output pairs of the system. In this case, every solution trajectory 
for the model (x, u, y) yields a solution (u, y) for the input-output equations, and every solution 
to the input-output equations corresponds to at least one model solution. 

In the usual differential algebra approach the input-output equations are generated as part of a 
characteristic set [111 I17| , but in principle input-output equations could be generated by a range of 
methods, e.g. simply by ad hoc substitution and differentiation, using Grobner bases (provided we 
add sufficient derivatives to generate enough equations), or via differential Grobner bases. Indeed, 
we show below that the coefficients of the input-output equations in general are identifiable, and 
that these coefficients yield complete identifiability information for the model, regardless of how 
they are generated. This may allow us to simplify and speed up the identifiability analysis process, 
as we can more flexibly tailor our methods to the system at hand. 

Theorem 4.1. The parameters of a rational function ODE model x(t, u, p), y(x, u, p) are globally 
(respectively locally) structurally identifiable if and only if the map c{p) from the parameters to the 
coefficients of a set of input-output equations is injective (respectively, the fibers contain finitely 
many elements), regardless of how the input-output equations are generated. 

Proof. Let 'i'(u, y, p) be the monic input-output equations of our system, with each %(u, y, p) := 
Ci(p)i/)i(u, y) = 0, where the Vi are distinct monomials in the inputs, outputs, and their 
derivatives, with the coefficient of the highest ranking term of each ipi equal to one. We refer 
to the complete set of coefficients of * as c(p). 

Let ^cs be a set of input output equations generated using the characteristic set approach of 
OUT]. Then if p and p* both yield the same trajectory y, the characteristic set approach tells 
us that the coefficients of ^cs satisfy ccs(p) = Ccs(p*)- Additionally, the results of Diop [H] 
and Saccomani [5] (with some algebraic manipulation) give us that the differential ideal generated 
by ^cs is prime, so that it coincides with its radical. As ^ and ^cs are both sets of input 
output equations, they have the same input-output solution pairs (i.e. the same variety). Then ^ 
is contained in the differential ideal generated by *cs (by Seidenberg's differential nullstellensatz, 
or by Hilbert's nullstellensatz after appending sufficiently many derivatives of Vl/cs)- Then as p 
and p* generate the same coefficients for ^cs, which in turn generates we know p and p* 
must both generate the same coefficients for yp, i.e. c(p) — c(p*). Then any solutions to the 
map $ : p — > y are also solutions to c(p), i.e. we have the forward implication that if c(p) is 
injective (respectively has finite fibers), then the model is globally (locally) identifiable, i.e. the 
map <& : p — >■ y is injective (respectively has finite fibers). 

To prove the converse, we need to show that if c(p) is not injective then neither is $, and if 
c(p) has infinitely many solutions, so does $(p). One way to show this is to show that if p and 
p* satisfy c(p) = c(p*), then p and p* also satisfy <i>(p) = <&(p*) (in other words, that p and p* 
both yield the same input-output pair, and so are indistinguishable). However, this follows from 
the definition of c(p), as if c(p) = c(p*), then by definition p and p* yield the same trajectory. 



5 



This means that solutions to the model for both p and p* yield the same input-output pair, i.e. 
$(p) = <f>(p*). □ 



This confirms that input-output equations contain complete identifiability information, which 
seems natural, as all the information obtainable from the known inputs and outputs (data) are 
present in the input-output equations. Next, we show that both Mansfield differential Grobner 
bases and Grobner bases (with sufficiently many derivatives of the model equations added in order 
to eliminate all non-observed variables) will yield a set of input-output equations. For this, we 
use the elimination theorem |15j and differential elimination theorem of Mansfield ^14 . This may 
restrict the type of ordering used, as Mansfield proved the differential elimination theorem using 
lexicographic order while Cox considers a braoder range of possible orderings for the non-differential 
version (see [TS] for more details). However, this is not a problem for identifiability analysis as 
lexicographic ordering is typically used [TJ [S] . 

Theorem 4.2. Both Mansfield differential Grobner bases and Grobner bases (with suffiently 
many derivatives of the model equations added) of a rational function ODE model x(i, u, p) = 
/(x, u, p), y(x, u, p) — g{x,t,p) contain a complete set of input output equations. 

Proof. Let us first show this for differential Grobner bases. Let / = x(t, u, p) — /(x, t,u,p), 
and g = y(x,u, p) — g{x,t,p), so that the model is given by / = 0, 5 = 0. Let mgb{f,g) be 
a Mansfield differential Grobner basis of the model (clearing denominators of model as needed), 
taking a lexicographic ordering with u<u<---<y<y<--- <x<x<--- (as given in [H]). 
Then {mgb} = {f,g}, so V{{mgb}) — V{{f,g}) and we have that solutions to mgb are precisely 
the allowed trajectories of the model. We also know that the characteristic set char{f,g) C 
{f,g} = {mgb}. We note that the characteristic set has been shown to contain a set of input 
output equations [SHHIIIT], so that {f,g} n R{y, u} 7^ 0. Then by the differential elimination 
theorem of Mansfield [13], we have that mgb D R{y, u} is a Mansfield Grobner basis for the ideal 
{/,5}nR{y,u}, i.e. 

{mgb n R{y, u}} ^ {/, g} n R{y, u}. 

Moreover, because solutions of {/, g} n R{y, u} are precisely the allowed input-output trajectories, 
and mgbn R{y,u} generates this ideal, then V{mgbr] R{y, u}) is precisely all input-output tra- 
jectories. Thus mgb contains a complete set of input output equations. A very similar argument 
shows that the same is true for algebraic Grobner bases, provided we begin by taking sufficiently 
many derivatives of the model equations to include all the differential variables needed to generate 
the characteristic set (and use the usual Grobner basis elimination theorem [15j rather than the 
differential version) . The algebraic Grobner basis result has also been shown independently in [T3] , 
with more details given there on the number of derivatives needed to generate the input output 
equations. □ 
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Figure 1: Linear 2-compartment model. 

5 Examples 

Next, we give some example applications of these results, and compare three of the methods for gen- 
erating the input-output equations discussed here — characteristic sets, differential Grobner bases, 
Grobner bases, and ad hoc substitution. All examples were implemented in Mathematica Version 8 
[22] . It must be emphasized that the speed of any of these methods depends heavily on the details 
of its implementation, so that results of speed tests such as these may vary widely. These examples 
are only intended to demonstrate that using different approaches can in some cases significantly 
improve the speed of the identifiability computations. 



Example 1: Linear 2-compartment Model. To begin, we introduce the linear 2-compartmcnt 
model, which is frequently used in pharmacokinetics, and has been shown previously by several 
methods to be unidentifiable [51 [TT] : 

Xl = u{t) + ki2X2 - {koi + k2i)xi 

±2 = k2iXi - (fco2 + ki2)x2 (2) 

y = xi/V 

where xi represents the mass of a substance in the blood (e.g. a hormone or drug), and X2 
represents the mass of the substance in the tissue. The drug exchanges between blood and tissues, 
and is degraded/lost in both compartments, at the rates given by the fc^j's above. The function 
u{t) represents a known input of the drug into the blood. The model output y = xi/V is the blood 
concentration of the drug, where V is the blood volume. The /c^j's and V are unknown parameters 
to be estimated. 

To examine the identifiabilty of ([2]), we start by generating a set of input-output equations. 



By Theorem 4.1 we need not use characteristic sets, so we will instead use ad hoc substitution to 



generate the input-output equations. We start by replacing xi with xi — yV , to give: 

yV = u{t) + ki2X2 - {koi + k2i)yV 

X2 = k2iyV - (fco2 + fci2)x2. 
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Next, we solve the first equation for X2, to yield X2 = ^W+''oiVy+k2iVy+Vv ^ ^j^jj,]-^ plug into 
the second equation (differentiating to give X2) to yield: 

ki2k2iVy + (fco2 + fci2) {u{t) - V{koi + k2i)y + Vy) + u{t) - V^(fcoi + k2i)y + Vy 

Collecting terms and making the polynomial monic (by dividing by the coefficient of the leading 
term y) yields the input output equation: 

- ^^^^ m(^) - + (fcoi^o2 + fcoifci2 + ^02^21)2/ + {km + fco2 + fci2 + k2i)y + y. (3) 

We note that this is the same input-output equation that would be achieved via a characteristic set, 
but requires fewer steps than the characteristic set algorithm (in part because the characteristic 
set maintains a set of three equations, each with different leaders, which must all be reduced with 
respect to one another, rather than focusing on just the input-output equation). However, this 
example is simple enough that both approaches take similar amounts of computational time (the 
substitution approach took slightly less time than the characteristic set, but the improvement was 
basically insignificant). 

In calculating the input-output equation, the differential Grobner basis approach took 0.022579 
seconds of CPU time in Mathematica, the substitution approach took 0.028624 seconds, and the 
characteristic set approach took 0.029163 seconds. Using a standard Grobner basis to calculate 
the input-output equations took significantly longer, at 0.339387 seconds of CPU time. 

By Theorem 4.1 the coefficients to ([s]) are identifiable. Then to test the identifiability of the 



individual parameters (fcoi, fco2, ^12, ^21, y)^ we must test injectivity of the map c(p). Thus, suppose 
we have an alternative set of parameters (01,02,03,04,05) which also yield the same output. As 
the coefficients for the input output equation are identifiable, we have that: 



^02 



ki2 02 + 03 



V 05 
1 _ 1 

koiko2 + fcoifci2 + ^02^1 = aia2 + aios + ^204 
koi + ko2 + ki2 + k2i = oi + 02 + 03 -f 04 

Solving for (fcoi, fco2, ^12, fei, V) following the algorithm in [llj reveals that the model is unidenti- 
fiable, with identifiable combinations fcoi + ^21, ^02 + ^12, ^12^21 and one identifiable parameter, 
V (as also shown in |11|). 



Example 2: Nonlinear 2-compartment Model. A common variant of the two compartment 
model is the following nonlinear version: 

il = u{t) + ki2X2 - ( -J^^ 1- ^21 ) Xl 

X2 = fc2lXi - (fco2 + fcl2)a;2 

y = xilV. 
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The input-output equation for this model is given by: 



V{ko2 + ki2) 2 2(fco2 + ki2) ko2 + ki2 Vy^ii 2yu ii ^02^21^^?/^ 



K'^ K V K"^ K V ' 

fc02V2+A:i2F2 + fc2lT/2\ 2. , ('^k^2k2lV ^k^2yVma. + k^2Vy^^a^\ 2 , '^^ VV 

P'^^^K P 

yy+ K02K21 H 7; 2/+ K02 + kx2 + ^21 + -7; — y+y 



^ . 1 uz, ^ . , UZ, ^ 

Solving for the parameters using the coefficients shows that the model is globally structurally 
identifiable. The timings for the calculation of the input-output equation for each of the differ- 
ent methods was: 0.013182 seconds CPU time in Mathematica for the differential Grobner basis 
approach, 0.164781 seconds for the substitution approach, 0.214768 seconds CPU time for the 
characteristic set approach, and 3.06737 seconds CPU time using Grobner bases. 



6 Conclusions 

In this work, we have extended the characteristic set-based differential algebra approach to struc- 
tural identifiability [51[TT1[T7] to incorporate more general differential polynomial solution methods. 
We show that the input output equations for a rational function ODE model, regardless of how 
they are generated, contain complete identifiability information for the model. In particular, we 
also show that Mansfield differential Grobner bases and algebraic Grobner bases can both be used 
to generate a set of input-output equations. These results may help to address one of the major 
criticisms of the differential algebra approach to identifiability — the computational complexity of 
the calculations required is often prohibitive for more complicated models [T]. 

For example, this is the case for the SIWR model examined in [23j, where the characteristic 
set approach was computationally infeasible, but using ad hoc substitution, the identifiability of 
a disease model used for cholera was able to be established. Indeed, in both examples here, the 
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substitution approach performed quite well (second fastest), suggesting that incorporating this into 
a more rigorous algorithm might improve computation time for identifiability. The substitution 
approach used here has the disadvantage that it is non-algorithmic (and so not general for other 
models), but it does allow one to fine tune calculations and speed up computation for a particular 
model. 

In the two examples given here, the differential Grobner basis approach gave the shortest com- 
putation time. The improvement in speed of differential Grobner bases over the algebraic Grobner 
basis approach may be due to the fact that algebraic Grobner bases do not take explicit advan- 
tage of the differential structure of the ring, making them less computationally efficient than the 
differential Grobner basis. However, as noted above, it must be emphasized that the details of the 
implementation for the Grobner basis algorithm may change these results significantly. Nonethe- 
less, this suggests that efficient algorithms for Mansfield differential Grobner basis calculation along 
the same line as those for algebraic Grobner bases (e.g. the Grobner walk [19]) might be a viable 
approach to improving the speed of identifiability calculations in practice. 

Acknowledgments 

This work was supported by National Science Foundation Award 0635561. Many thanks to Joseph 
Tien, Suzanne Robertson, Jeff Dunworth, and Tony Nance (Ohio State University) for their helpful 
comments. 

References 

[1] G. BcUu, M. P. Saccomani, S. Audoly, and L. D'Angio. Daisy: a new software tool to test global 
identifiability of biological and physiological systems. Comput Methods Programs Biomed, 
88(1):52-61, 2007. 

[2] Michael J. Chappell and Roger N. Gunn. A procedure for generating locally identifiable repa- 
rameterisations of unidentifiable non- linear systems by the similarity transformation approach. 
Mathematical Biosciences, 148(1):21-41, 1998. 

[3] C. Cobelli and J. J. DiStefano. Parameter and structural identifiability concepts and ambigu- 
ities: a critical review and analysis. American Journal of Physiology - Regulatory, Integrative 
and Comparative Physiology, 239(1) :R7-R24, 1980. 

[4] H. Pohjanpalo. System identifiability based on the power series expansion of the solution. 
Mathematical Biosciences, 41(l-2):21-33, 1978. 

[5] Maria Pia Saccomani, Stefania Audoly, and Leontina D'Angi. Parameter identifiability of 
nonlinear systems: the role of initial conditions. Automatica, 39(4):619-632, 2003. 

[6] R. Bellman and K. J. Astrom. On structural identifiability. Mathematical Biosciences, 7(3- 
4):329-339, 1970. 

[7] J. J. Distefano. Complete parameter bounds and quasiidentifiability conditions for a class of 
unidentifiable linear systems. Mathematical Biosciences, 65(l):51-68, 1983. 



10 



[8] L. Ljung and T. Glad. On global identifiability for arbitrary model parameterization. Auto- 
matica, 30(2):265-276, 1994. 

[9] N. Meshkat, M. Eisenberg, and J. J. Distefano. An algorithm for finding globally identifi- 
able parameter combinations of nonlinear ode models using Groebner bases. Math Biosci, 
222(2) :61- 72, 2009. 

[10] J. F. Ritt. Differential Algebra. New York: Academic Press, 1950. 

[11] S. Audoly, G. Bellu, L. D'Angio, M. P. Saccomani, and C. Cobelli. Global identifiability of 
nonlinear models of biological systems. IEEE Trans Biomed Eng, 48(1) :55 65, 2001. 

[12] Gabriella Margaria, Eva Riccomagno, Michael J. Chappell, and Henry P. Wynn. Differential 
algebra methods for the study of the structural identifiability of rational function state-space 
models in the biosciences. Mathematical Biosciences, 174(1):1 - 26, 2001. 

[13] N. Meshkat, C. Anderson, and J. J. DiStcfano. III. Alternative to Ritt's Pscudodivision 
for finding the input-output equations in algebraic structural identifiability analysis. ArXiv 
e-prints, March 2012. 

[14] Elizabeth Mansfield. Differential Grohner Bases. PhD thesis, University of Sydney, 1991. 

[15] D. O'Shea D. Cox, J. Little. Ideals, Varieties, and Algorithms: An Introduction to Computa- 
tional Algebraic Geometry and Gommutative Algebra. Springer, 1996. 

[16] Christian Aistlcitner. Relations between grobner bases, differential grobner bases, and differ- 
ential characteristic sets. Thesis, Institut fiir Symbolisches Rechnen, 2010. 

[17] F. OUivier. Le problem de I'identifiabilite structurelle globale: etude theoretique methodes 
effectives et bornes de complexite. PhD Thesis, Ecole Polytechnique, 1990. 

[18] Jean-Charles Faugerc. A new efficient algorithm for computing groebner bases (f4). Journal 
of Pure and Applied Algebra, 139(1-3):61 - 88, 1999. 

[19] S. CoUart, M. Kalkbrener, and D. Mall. The grbner walk. Technical report, J. Symbolic 
Computation, 1993. 

[20] G. Carra'Ferro. Groebner bases and differential algebra. In Proceedings of the 5th international 
conference, AAEGG-5 on Applied Algebra, Algebraic Algorithms and Error- Gorrecting Godes, 
pages 129-140, New York, NY, USA, 1987. Springer- Verlag New York, Inc. 

[21] Sette Diop. Differential-algebraic decision methods and some applications to system theory. 
Theor. Gomput. Sci, 98(1):137-161, May 1992. 

[22] Wolfram Research, Inc. Mathematica edition: Version 8.0. Wolfram Research, Inc., Gham- 
paign, Illinois, 2010. 

[23] M. Eisenberg, S. Robertson, and J. Tien. Identifiability and estimation of multiple transmis- 
sion pathways in cholera and waterborne disease. Journal of Theoretical Biology, in press. 



11 



